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REVISION NOTICE 

This description replaces previous descriptions of "Matrix Inversion 4," 
program Dl- 139.0. Program references have been changed to current 
designations. 



FUNCTION 

"Matrix Inversion 4" enables the source program to replace the elements 
of a square matrix with the elements of its inverse. Matrices of any 
rank greater than one will be accommodated. The routine is entered and 
left in machine language, but it uses Floating Point Interpretive 
System 4, program Hl-24.3 for all calculations. 



INPUT 

The elements of a square matrix are stored in consecutive locations on 
the memory drum beginning in Mo. The elements must be in double 
precision floating point format. 



OUTPUT 

The elements of the inverse matrix are stored in consecutive locations, 
beginning in Mo, in double precision floating point format. 
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PRINTED IN U.S.A. 



MATRIX INVERSION 4 



CALLING SEQUENCE 
Location 



Order 

E 
R 
U 
Z 
Z 
Z 
etc. 



Address 

0000 

Lo 

Lo 

Lo of source program 

Mo 

N 



The E0000 order in (a - 1) is required only if the previous instruc- 
tions are interpreted by the source program. N, in (a + 4), is the 
rank of the matrix. 



TIME 



Approximately 2.2. N^ milliseconds are required. 



STORAGE 

2 tracks and 13 sectors (141 words) are required in memory for in- 
structions and constants. No temporary storage is used. Although 
only N sectors are required for the matrix elements, this routine 
requires N^ + N sectors beginning in Mo. 



NOTES 

Matrices of rank greater than 30 can not be inverted due to machine 
storage limitations. 
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FUNCTION 

To replace the elements of a square matrix with the elenx-n-a o* its 
inverse. Matrices of any rank greater than one vill be accomodated, 2he routine 
is entered and left in machine language, but it uses the Double Pre- J slot. 
Floating Point (DPFP) program ( ), for all calculations. 

IMHJT 

The elements of a square matrix in consecutive locations on the dnm. 
beginning in Mo. The elements most be in Double Precision Floating Point farm. 

OUTPUT 

^ . J^l elements of the inverse matrix in consecutive locations, beginning in 
mo, in DtrF form. 



CALLING SBtyJEBHCS 






Location 


Order 

S 


Address 


a - 1 


0000 


a 


R 


Lo 


a + 1 


U 


Lo 


a + 2 


z 


Lo of I 


a ♦ 3 


z 


Mo 


a + k 


z. 


n 


a + 5 


etc. 





Die EOOOO order in a - 1 is required only if the previous instructions 
are interpreted by the DPFP program. N, in a + k, is the rank of the matrix. 

TIME 

Approximately 2.2 n seconds. 

STORAGE 

Two tracks and 13 sectors are required for instructions and constants. 
Mo temporary storage is used. Although only n^ sectors are required for the matrix 
elements, the routine requires n x ♦ n sectors beginning in Mo. 

NOTES 

Matrices of rank greater than 50 can not be inverted due to machine storage 
limitations • 
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